Nonlinearity-induced conformational instability and dynamics of biopolymers 

Serge F. Mingaleev^'^"'^, Yuri B. Gaididei^''^, Peter L. Christiansen^, and Yuri S. Kivshar^ 
^ Nonlinear Physics Group, Research School of Physical Sciences and Engineering, 
Australian National University, Canberra ACT 0200, Australia 
^ Department of Mathematical Modeling, The Technical University of Denmark, DK-2800 Lynghy, Denmark 
^ Bogolyuhov Institute for Theoretical Physics, 14-B Metrologichna St., 03143 Kiev, Ukraine 

We propose a simple phenomenological model for describing the conformational dynamics of 
biopolymers via the nonlinearity-induced buckling and collapse (i.e. coiling up) instabilities. We 
describe the buckling instability analytically, and then demonstrate the role of both instabilities in 
the folding of biopolymers through the numerical simulations of the three-dimensional dynamics of 
a long semiflexible chain in the aqueous environment. 
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Conformational flexibility is a fundamental property of 
polymers which differentiates them from small molecules 
and gives rise to their remarkable properties. However, 
even the properties of polymers seem meagre when com- 
pared to the functionality of biopolymers. A distinctive 
feature of biopolymers is that, on the one hand, they are 
heterogeneous and, on the other, their elementary sub- 
units have a complex structure and can carry long-lived 
nonlinear excitations . While the significance of the 
heterogeneity gained wide recognition as being of funda- 
mental importance for protein folding the function 
of nonlinear excitations remains obscured, in spite of a 
large body of theoretical research [^-Q. 

Here we suggest a new role of nonlinear excitations 
in one of the most important functionalities of biopoly- 
mers, their conformational dynamics. For the first time 
to our knowledge, we demonstrate that nonlinear exci- 
tations can facilitate bending or folding of a semiflexible 
molecular chain via the buckling or collapse instability, 
and thus provide a possible scenario for the conforma- 
tional dynamics of biopolymers. We consider a semiflex- 
ible chain that models a biopolymer and assume that its 
mechanical degrees of freedom are coupled to the internal 
degrees of freedom such as amide-I vibrations in proteins 
base-pair vibrations in DNA 1^,^, or polarons in 
proteins |^,^ . This coupling provides the energy needed 
to overcome the potential barriers for folding, thus avoid- 
ing the need for stochastic thermal activation. We reveal 
that such a simple model describes the conformational 
dynamics as a deterministic controlled process: when the 
amplitude of the internal excitation becomes large, the 
chain displays the buckling or collapse instability that 
subsequently initiates the chain folding. We believe that 
such simple physics can also account for the kinetics of 
conformational phase transitions of semiflexible polymers 
in solutions 

The importance of nonlinear excitations has already 
been emphasized in the modeling of biopolymers and con- 
jugated polymers. In particular, it was suggested that 
solitons may provide a possible physical mechanism for 
the energy (or charge) transport and storage in proteins 
and conducting polymers |^ . The concept of nonlin- 



ear localized modes was also employed to explain some 
specific features of the DNA dynamics [^|j4j. However, 
most of those studies were dealing with nonlinear prop- 
erties of straight molecular chains and never discussed 
changes in their geometrical structure. The study of non- 
linear properties of curved chains has been initiated only 
recently [ p^ and, in particular, it was shown that the 
bending in a curved chain manifests itself as an effec- 
tive trap for nonlinear localized modes. In this Letter, 
we make an important step forward and, taking into ac- 
count flexibility of biopolymer chains, demonstrate that 
nonlinear excitations may act as drivers giving impetus 
to conformational dynamics of biopolymers. 

Let us consider a simple phenomenological model of 
a biopolymer consisting of particles with mass M, in- 
dexed by n, and located at the points r„ = {xn, Vn, Zn}. 
We assume that each particle represents a complex sub- 
unit of the biopolymer which, additionally to its position 
r„ (t) , carries an internal excitation which can be charac- 
terized, in some approximation |^-||,^, by the complex 
amplitude, ?A„(t). Such an internal mode can represent, 
for instance, a polaron state (an excess electron accompa- 
nied by phonons) or vibrational state (amide-I vibrations 
in proteins or base-pair vibrations in DNA). The Hamil- 
tonian of such a chain is written as iJ = T+U + V, where 
T = iM/2)J2nid-rn/dty is the kinetic energy, C/(r„) is 
the potential energy of inter-particle interactions, and 

n I m^n J 

is the energy of the internal excitation. Here x character- 
izes on-site self-trapping nonlinearity of the internal ex- 
citation, and Jnm are the excitation transfer coefficients 
that depend on the distance in the embedding space be- 
tween the particles n and to: Jnm = J{\'''n ~ ''ml) — 
(e" — 1) exp (— Q;|r„ — r,„|) From the Hamiltonian, 
we obtain the equations of motion 

k m^k 
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where we allowed additionally for a viscous damping v of 
the aqueous environment. We assume that the number 
of excitations, J2n \^n\^ = 1, is conserved by Eq. (j^). 
This approximation is fully justified for polaron states 
and should be fairly good for vibrational states, taking 
into account that in the case under study the lifetime of 
amide I vibrations is greatly inceased as a result of their 
self-trapping locahzation [|]J|]. 

We consider the potential energy of inter-particle inter- 
actions as a sum, [/(/•„) = Us+Ub+Ur, of the stretching 
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inextensible (u oo), so that the distance between the 
neighboring sites of the chain remains constant. In this 
case, a spatial conformation of the chain is completely 
defined by the bending angles On- 
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and (to describe correctly coiling up of the chain) the 
energy of short-range repulsive interactions between par- 
ticles (considered as elastic balls of the diameter d): 
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for |r„ — rra\ < d, and Ur — 0, otherwise. Here a is the 
equilibrium lattice spacing (we take a = 1), cr is the elas- 
tic modulus of the stretching rigidity of the chain, f?„ is 
the angle between neighboring bond vectors, (r„ — r„_i) 
and (r„+i — r„), that meet at the n-th particle [see Fig. 
|l|(a)], ^max is a maximum bending angle, and k is the 
elastic modulus of the bending rigidity of the chain. 

First of all, we recall the properties of nonlinear local- 
ized modes in a curved chain, when the chain geometry is 
"frozen" . Such an analysis has been recently carried out 
in Ref. , where it was shown that the chain curvature 
affects strongly the properties of nonlinear modes, creat- 
ing an effective double- well potential around the bending 
point, which is responsible for the symmetry-breaking ef- 
fect: as the nonlinearity x increases above some thresh- 
old, a symmetric stationary mode becomes unstable, and 
it transforms into an energetically more favorable asym- 
metric stationary mode. In that study, it was also found 
that the energy of a nonlinear mode decreases as a square- 
law of the curvature when the curvature is small. Thus, 
since the bending energy of the chain increases also as the 
square of the curvature, one can expect that there should 
exist a converse influence of the internal nonlinear modes 
on the conformational dynamics of the chain. Below, we 
confirm this conjecture by straightforward numerical cal- 
culations and find that not only the chain tends to bend, 
for small enough k, but it also collapses (i.e. coils up), 
for smaller k. 

As the first step of our study in this Letter, we analyze 
the stationary ground-state configurations of the chain 
described by the model introduced above. We assume 
that the chain is infinitely long, planar (z„ = 0), and 



FIG. 1. (a) Schematic representation of a buckled chain. 
The distance between the sites near the bend becomes smaller 
(-R2 < -Ri) lowering the coupling energy of a nonlinear mode 
that works against the increased bending energy and lead to 
a buckling instability, (b) Ground state spatial configurations 
of a semiflexible chain (a — 2, 6'max = 7r/3, and k — 0.125) 
at different values of x: 2.13 (i), 2.26 (ii), 2.44 (iii), 2.88 (iv), 
3.00 (iv), 3.13 (iii), 3.18 (ii), and 3.25 (i). As the nonlinearity 
X increases, the chain bending initially increases, reaches its 
maximum at x ~ 2.96, and then decreases. 
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FIG. 2. Phase diagram of the ground states of a semiflexi- 
ble biopolymer chain calculated numerically (solid curves with 
squares) and analytically from Eqs. (^-(|^ (dashed curve), at 
a — 2, Smax = 7r/3, and u = 00. 

We find the ground-state configurations of the chain 
looking for the stationary solutions of Eq. (|^) in the 
form ^/jn{t) = e^^^cfin and drn/dt = 0. The corresponding 
system of nonlinear algebraic equations has been solved 
by the Newton-Raphson iteration scheme. The results 
of these calculations are presented in Fig. ^ The most 
interesting feature of this phase diagram is the existence 
of three distinct types of the ground states. 

A stiff chain (with n > 0.183) is always straight and, 
regardless of the value of x, it cannot be affected by 
the presence of a nonlinear excitation. However, for 
K < 0.183, there appears a finite domain of x, for which 
a straight chain becomes unstable in the presence of a 
nonlinear mode. Inside of this domain, the chain in 



2 



its ground state curves symmetrically around the non- 
linear excitation, with some finite overall buckling angle 
/3 = On [see Fig. |](b)]. In Fig. ^, we plot the depen- 
dence of the buckling angle f3 vs. x for different values 
of K. As is seen, for 0.124 < k < 0.183 the buckling angle 
f3 is always finite and reaches its maximum at x ~ 3. 

When the bending rigidity of the chain decreases fur- 
ther, i.e. for K < 0.124, a new type of the ground state, 
a collapsed chain, emerges (see Figs. ^ and ||). In this 
case, the distance between the complementary particles 
of the chain [see Fig. |l](b)] decreases to an extent that 
the growing interaction between the tails of the nonlinear 
excitation turns buckling into collapse causing the coil- 
ing up of the chain. The conformation of the chain in 
the collapsed state (which is usually a three-dimensional 
compact globule, see Fig. H) is sensitive to the peculiar- 
ities of the potentials (121)-^) and initial conditions. 
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FIG. 3. Buckling angle (3 vs. the nonlinearity x for a = 2, 
Omax = 7i"/3, a = oo, and varying k. 

To gain a better insight into the physical mechanism 
of the buckling instability, we have developed a simple 
variational approach. We assume that the chain bending 
occurs only at the center of the excitation (n = 0). In this 
case, the chain is composed of two semi-infinite lines with 
the overall buckling angle (3 = Bq between them. Obvi- 
ously, such a geometry cannot describe the chain collapse 
because /3 is restricted by the maximum angle ^maxj but 
it should be valid as long as we are not interested in the 
collapsed phase but study only a threshold between the 
straight and buckled chain phases. 

In such a simplified model, the ground state of the 
chain is determined by the dependence of the system en- 
ergy on (3. We use a variational ansatz for the localized 
mode in the form 0„ = \/tanh7 exp(— 7 |n|), which sat- 
isfies the constraint X^n ~ f i ^^'^ minimize H with 
respect to P and 7. Thus, we can show that the straight 
configuration of the chain (/3 = 0) corresponds to the en- 
ergy maximum (i.e., the straight chain is unstable) for 
K < Kcr, where the critical value of the bending stiffness 
is given by the following analytical expression 

3 e'^+'T — 1 

Kcr = XQ(e° - 1) 3(^0+7 _ 1)3 ^^^^^ ' f^"^) 

and 7 is determined from the equations 
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In Fig. 1^, we compare the threshold curve that separates 
the phases of the straight and buckled conformations pre- 
dicted by the analytical formula (Q)-(^, with the same 
dependence calculated numerically. One can see that the 
variational method provides reasonably fair analytical re- 
sults for the phase diagram. 

An important question is how the buckling and col- 
lapse instabilities manifest themselves in the dynamics 
of semiflexible biopolymer chains. To clarify this issue, 
we consider a simple (but physically important) case of 
an initially straight chain of a finite length whose left end 
is free whereas the right end is fixed. First, we excite a 
single particle at the left end of the chain, and then keep 
track of how this excitation propagates along the chain 
and how it modifies the chain geometry. For small val- 
ues of X, such a single-particle excitation disperses and 
spreads out, whereas for very large values of x it does not 
propagate and remains completely trapped at the end of 
the chain. However, for the intermediate values of x (in 
our model, for 2.5 < X < 3.5), such a single-particle 
drive applied at the end of the chain generates a moving 
localized mode. It is remarkable that this interval of non- 
linearity parameter x matches closely the region in which 
the bending capability of the nonlinear mode reaches its 
maximum. 

For a stiff biopolymer chain (i.e., for k > 0.183), the 
propagation of the generated nonlinear mode does not 
produce any conformational change of the chain. How- 
ever, when K decreases, the nonlinear mode generates 
a local chain bending. In this case the nonlinear mode 
propagates along the chain being accompanied by a lo- 
cal bending; the amplitude of such a bending decreases 
as the viscous damping v increases. Thus, the buckling 
instability we have described above appears to be physi- 
cally visible, and it is responsible, in particular, for large- 
amplitude localized bending waves. 

Let us consider now a flexible biopolymer chain with 
even smaller k, for which the nonlinear mode should pro- 
duce a collapse instability (see Fig. ||). In this case, the 
evolution of the chain clearly demonstrates a collapse dy- 
namics (see Fig. ^), and the chain folds into a compact 
coil with several loops consisting of 7-10 particles each 
p3[ . Importantly, even for large values of the damping 
parameter v the collapse of the chain is only delayed, but 
not prevented. For instance, for the parameters used in 
Fig. U the collapse-on time tc grows almost linearly with 
ly: tc « 87+ 1290:/+283j/2 for < 1. Besides, increasing 
of ly leads to decreasing of the distance between the place 
of the collapse nucleation and the left end of the chain. 

It should be mentioned that the instability we discuss 
in this Letter remains latent in a straight infinitely long 
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chain, because the bending of such a chain would require 
an infinite energy. However, as we have demonstrated 
above, the instabihty manifests itself as soon as we con- 
sider more realistic cases and take into account a finite 
length of the chain. In this case, a nonlinear localized 
mode, once generated at the end of the chain, propa- 
gates along the chain accompanied by a local chain bend- 
ing, for intermediate values of the bending rigidity k, and 
it causes the chain folding, for smaller k. 
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FIG. 4. Dynamics of a semiflexible biopolymer chain due to 
the collapse instability produced by a moving nonlinear mode 
generated from the single-particle excitation at the left end of 
the chain. The chain consists of 30 particles, with the right 
end fixed. Marked particles indicate location of the moving 
excitation (for these particles, lipnf > 0.25). The parameters 
are: x = 3.2, k = 0.06, u = 0.3, a = 10^ a = 2, 0max = 7r/3, 
M = 0.5, S = 10^ and d = 0.6. 

In conclusion, we have predicted analytically and 
demonstrated numerically a new role played by nonlinear 
excitations in the conformational dynamics of biopoly- 
mers. Taking into account the coupling between the in- 
ternal and mechanical degrees of freedom of a semiflexible 
biopolymer chain, we have found that it may produce the 
buckling and collapse instabilities of an initially straight 
chain. In the case of collapse instability, the viscous 
damping of the aqueous environment only slows down 
the folding (coiling up) of the chain, but does not stop 
it even for relatively large damping (z^ ~ 1). It should 
be stressed that these instabilities are most pronounced 
for intermediate values of the self-trapping nonlinearity, 
and they vanish completely in the linear limit (i.e., for 
X = 0) . We find that this effect is only weakly affected by 
the peculiarities of the interaction potentials , and 

thus it should be generic for different models of nonlinear 
semiflexible chains. 

We believe that the nonlinearity-induced buckling and 
collapse instabilities may be of a considerable importance 
in the conformational dynamics of the globular macro- 



molecules with a multitude of quasi-isoenergetic con- 
formations 1^, where nonlinear modes excited even far 
from the chain ends should cause conformational transi- 
tions between different types of the globular states of the 
same or similar topology. Presumably, another related 
problem is the energy transduction in molecular motors, 
for which it was recently suggested that the free energy 
from ATP hydrolysis is initially converted into the hinge- 
bending motion of the motor sub-unit p4[ . 
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